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Abstract 

We study the critical behavior of the systems dominated by Coulom- 
bic interaction. For this purpose we used the method of collective 
variables with a reference system. Starting from the Hamiltonian of 
the restricted primitive model (RPM), the simplest model of ionic 
fluids, we obtain the functional of the grand partition function given 
in terms of the two collective variables: the collective variables pk 
and Ck describing fluctuations of the total number density and charge 
density, respectively. As the result of integration over the variables 
Ck, the microscopic based effective Hamiltonian of the RPM at the 
vicinity of its gas-liquid critical point was constructed. The coeffi- 
cients of the effective Hamiltonian, describing the density fluctuations 
nearby the gas-liquid critical point, are analyzed. It is shown that, 
in spite of the long-range character of the Coulombic potential, the 
effective interactions appeared at this level of the description have 
the short-range character. As a result, the effective Hamiltonian ob- 
tained for the RPM in the vicinity of the critical point has the form 
of the Ginzburg-Landau- Wilson Hamiltonian of an Ising-like model in 
a magnetic field. This confirms the fact that the critical behavior of 
the RPM nearby the gas-liquid critical point belongs to the universal 
class of a 3D Ising model. 
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1 Introduction 



Experimental investigations of the critical properties of electrolyte solutions 
displayed three different types of behavior: Ising-like and mean-field behavior 
as well as a crossover between the two PQ-|E]- In order to interpret the results, 
ionic systems were classified in solvophobic and Coulombic. In solvophobic 
systems, Coulombic forces are not supposed to play a major role; the critical 
behavior is that common for fluids and fluid mixtures, i.e. Ising-like. By con- 
trast, in Coulombic systems the phase separation is driven by electrostatic in- 
teractions. In recent years the critical behavior of the Coulombic systems has 
been a subject of active research. A theoretical model which demonstrates 
the phase separation driven exclusively by Coulombic forces is the restricted 
primitive model (RPM) In this model the ionic fluid is described as an 

equimolar-equisized mixture of opposite charged hard spheres immersed in a 
dielectric medium. Early studies established [Zj-0 that the model has a gas- 
liquid (GL) phase transition. A reasonable theoretical description of the GL 
critical point in the RPM was accomplished at a mean-field (MF) level using 
integral equation methods EI] and Debye-Hiickel theory JI] . However, 
a systematic theoretical investigation of the criticality of the RPM requires 
the microscopic based effective Ginzburg-Landau- Wilson (GLW) Hamilto- 
nian. Although some progress in this direction has been made ^2]-CE] the 
GLW functional of the continuous RPM explicitly related to the microscopic 
characteristics has not been derived yet. 

For the last decade, the GL critical point in the RPM has been much 
studied by computer simulation methods [Hj-[Sn]. While the theory has not 
provided the clear picture of thermodynamic behavior of the model in the 
critical region,very recent simulation studies have found strong evidence for 
Ising universal class [23 HH] ■ 

In this paper we address an issue of the criticality in the RPM using 
the functional representation of the grand partition function in terms of the 
collective variables (CVs). The transformation from the individual to the 
collective coordinates (the Fourier transforms of the density fluctuations) is 
carried out via the corresponding Jacobian. Within the random phase ap- 
proximation for the Jacobian one arrives at the Debye-Hiickel approximation 
for the free energy j3T|. First, the CVs were introduced for the description of 
the charged particle systems in the 1950s The further development 

of the CV method was related to the theory of phase transition of the second 
type proposed in |33]-|57j. 
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According to the basic idea of universality [HHj, different systems which 
can be described by similar effective GLW Hamiltonians (the Hamiltonians 
of the same symmetry) demonstrate the same critical behavior. Thus, the 
knowledge of the effective Hamiltonian of the system it is important in de- 
terming its critical behavior. Our purpose here is to derive, from the first 
principles, the effective GLW Hamiltonian of the RPM in the vicinity of its 
GL critical point. 

The layout of the paper is as follows. We give a functional representation 
of the grand partition function of the RPM in Section 2. In Section 3 we con- 
struct the effective Hamiltonian. On this basis, some conclusions concerning 
the critical behavior of the RPM nearby the GL critical point are made in 
Section 3. 



2 Functional representation of the grand 
partition function of the RPM 

The RPM consists of N = N + + AL hard spheres of diameter a with iV + 
carrying charges +q and iV_ (= N + ) charges —q, in a medium of dielectric 
constant D. The interaction potential of the RPM has the form 



U l5 (r) 



oo 

Dr 



if r < a 
if r > a 



Qi = ±q- 



(1) 



We start with the grand partition function for a two-component system 
(7,<5 =+,-): 



5 = £ £ n 

N+>0 N->0f=+,- 

where (dT) = Y\ y dT N dT 



AT-, 



iV 7 ! 



A'. 



(dT) exp 



drldrl 



7^ 



75 ij 



dr]f (7 = +, — ) is an element 
of the configurational space of the 7th species; z 1 is the fugacity of the 
7th species: z 7 = exp(/3/x ), \i — fj, 7 + (3~ l ln[(27rm 7 /3 -1 ) 3//2 //i 3 ], /x 7 is the 
chemical potential of the 7th species and \x is determined from the equation 



<91nH 

<9/3/i 7 



(N 7 ). 



For the RPM we have fi + = fi. 
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Now we present the interaction potential U^{r) as a sum of two terms: 



U lS {r) = W(r) + $ 75 (r), 

where Vv( r ) * s a potential of a short-range repulsion and 3> 7 a(r) is a long- 
range attractive part of the potential. We split the potential U 7 g(r) into 
short- and long-range parts using the Weeks- Chandler- Andersen partition 
As a result, we have 



(r) 



oo if r<(7 
if r > a ' 



if r < cr 



'/-.. '/,> 
Dcr 



This simple form for $ 7 5(r) inside the hard core changes the behavior of 
the Fourier transform for large k from usual Coulombic k~ 2 to k" 2, decay. As 
was shown [40J, this choice of $ 7 5(r) for r < a produces rapid convergence 
of the series of the perturbation theory for the free energy. The Fourier 
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transform of $ 7 5(r) = |p = $c( r ) has the form 

/3p*c{x) = 24/f^, (2) 

where (3* = jf^, (3 = r\ = fpcr 3 is fraction density, x = ka. 

Within the framework of the approach considered the interaction con- 
nected with a repulsion (potential xjj y s( r )) is described in the space of the 
Cartesian coordinates of the particles. We call the hard sphere system with 
the diameter a a reference system (RS). The interaction connected with an 
attraction (potential $ 7 <s(r) ) is considered in the CV space. The transfor- 
mation from the Cartesian coordinates to the CVs is performed by means of 
the transition Jacobian. 

Using the method of CVs developed for a two-component continuous sys- 
tem jHJ-112] we can rewrite the grand partition function of the RPM in the 
following form: 

S = S J (dp)(dc) exp f /fyiPo - — ^ c k c_ k j J(p, c). (3) 
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Here S is the grand partition function of the RS. pk and Ck are the CVs 
which describe total density fluctuation modes and charge density fluctuation 
modes, respectively: 

Pk = p k - iPk> c k = 4 - ic k , 

the indices c and s denote the real and imaginary parts of CVs pk and Ck- 
Each of p k (c k ) and p k (c k ) takes all the real values from — oo to +00, and 
(dp) and (dc) are volume elements of the CV phase space: 

/ / 

(dp) = dpoY\ dp^dp^, (dc) = dc Y\_ cfc k cfc k . 

k^O k^O 

Unknown parameter p 1 (pi = (pi i+ + p\-)/\/2) is determined from the 
equation 

ainH i= (A0 

d(3p 1 y/2' V ' 

J(p, c) is the Jacobian of the transition to the CVs averaged over the RS. 
For the RPM, J(p, c) is of the same form as that for the symmetrical binary 
fluid 



J(p,c)= f (du) (d-f) exp [i27r^(wkPk + 7kC k ) + ^^ J Di i " ) (^ ) 7)], (5) 

" 1, m >*i A ^.n 



71>1 2n>0 



n! z — ' 

ki...k„ 

Tkx--- 7k In ^k In+1 • • • Wk„4 1+ ...+k n , (6) 



and variable c<Jk (7k) is conjugate to CV pk (c k ). Index z n (i n = 0, 2, 4, . . . 2n) 
is used to indicate the number of variables 7k in the cumulant expansion 
®. Cumulants Wl ( n n) are expressed as linear combinations of the partial 
cumulants 9Jt 7l ... 7n and are presented for 71, . . . , 7 n = +, — and n < 4 in [HI] 
(see Appendix B in [57]). 
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In ©-© the cumulants M { n n) with i n = are connected with the nth 
structure factors of the RS [37] : 

SOT?) = (JV)5 B . (8) 

Cumulants with i n ^ can be expressed in terms of 9Jt^ (see also formulae 
(4.8) in 

= 15071^3 - 309Jli°2 4 + 1697li°2 5 . (9) 

In general, the dependence of 9Jtn" (ki, . . . , k n ) on wave vectors ki, . . . , k„ 
is complicated. Since we are interested in the critical behavior, the small- 
k expansion of the cumulant can be considered. Hereafter we shall replace 
by their values in long-wave length limit 0tt£ n) (0, . . . , 0). 
Structure factors S n (0, . . . , 0) with n > 2 can be obtained from S^O) by 
means of a chain of equations for correlation functions |43j . 
Let us present J(p,c) as 

J (p, c) = / (dw) (dj) exp [i27r ^(u; k p k + 7 k c k ) + 



-«2tt 



,2 



k 

J2Y, D n n) (^)l (10) 



n>3 i n >0 

where Dn {^il) has the following form (up to n — 4): 

DS»\u,-y) = ^ (^ 0) ^ kl ^ k2 ^ k3 +3^ 2) a; kl 7 k2 7k3)4 1+ k 2+ k3 + 

ki,k 2 ,k 3 

C-i27r) 4 ^ - (o) «cm( 2 ) 

2 , (™4 ^k^kaWkgO;^ + om\ u; kl ^ k2 7 k3 7 k4 + 



4! 

ki,...,k 4 

>(4) 



9^4 7k 1 7k 2 7k 3 7k 4 )^k 1 +...+k4- 

In (jlOj) the linear term is eliminated by the shift p k = p' k + 9Jl[ 5 k (the prime 
on p k is omitted for clarity). 
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According to ©-©, Wl { 2 \k) = ^S 2 (k) and 9Jlf = ^, where S 2 (k) 
is a two-particle structure factor of a one-component hard sphere system. In 
the Carnahan-Starling approximation we have for S 2 (k = 0): 

S 2 = -. 

1 + 4r/ + 4r/ 2 - 4rf + r/ 4 

3 The effective Hamiltonian of the RPM in 
the vicinity of the GL critical point 

Because 9Jl 2 °\ty * s the positive and smooth function in the region under 
consideration and DJl 2 2 ^ is equal to constant, we can integrate in (fTUj) over a; k 
and 7k using the Gaussian density measures as basic ones. This integration 
can be performed by several ways. For example, using the Euler equations 
we can determine (and 7^) which provide the maximum for the functional 
in the exponent in 1)100: 

* P— k * c -k /1 1 \ 

^k = ^THT + • • • > 7k = ^TrT + (H) 

-i27r^ 0) -i2vr^ 2) 

As a result, we can present H in the form: 

H = S C / (dp)(dc) exp ofVo ~ ^[ (4 0) PkP-k + 4 2) c k c_ k ^ + 

' k 

1 ( (°) 1 ( 2 ) Ax 

7^ 2^ P-kiP-k 2 P-k 3 + a 3 P-k!C-k2 C -k 3 y fki+ka+ks + 

ki,k2,k3 

1 ( (°) , (2) 

4J P-kiP-k 2 P-k 3 P-k 4 + «4 P-k!P-k 2 C-k3C-k 4 + 



ki,...,k4 

,( 4 ), 



a\ ; c_ kl c_k 2 c_ k3 c_ k4 J 5 kl+ ...+k 4 (12) 
and super index i n indicates the number of variables at oi'"' '. Here 

c= n^wii^exp(^)), (13) 

and the coefficients On , taking into account the first terms in (jll)) . can be 
written as 



(o)_, fl( o)_ 1 a (2)_ 1 + 2p c (*QM? ) 
ij -ppi, a 2 - -— ^y, a 2 , (11) 



Wl v 2 > ~ m 2 



Another way of integration in (fTU|) is outlined in Appendix. 

First, we restrict our attention to the Gaussian approximation which 
corresponds to neglecting the terms proportional to p 3 , pc 2 , etc. in the 
exponent of (|12j) . As was shown |33], this approximation yields only the 
boundary of stability with respect to fluctuations of the charge density. The 
fact that the RPM does not demonstrate the GL phase instability in this 
approximatuin is attributed to the effects of charge-charge correlations not 
present at this level. In order to obtain the GL spinodal curve we should 
take into consideration the terms of the order higher than the second one 

mam. 

Now we follow the programme proposed in [37J for a two-component 
fluid system. First, we separate the two types of variables: the essential 
variables (which include the variable connected with the order parameter) 
and the non-essential variables. Then, integrating over the non-essential 
variables with the Gaussian density measure, we construct the basic density 
measure (the GLW Hamiltonian) with respect to the essential variables. 

For the RPM in the vicinity of the GL critical point the variable p k 
(describing the fluctuations of the total number density) turns out to be the 
essential variable [12] ■ Thus, we can present (|T2|) as 

1 v - 1 



E C J (dp) exp a lPo + ^ a 2 p k p_ k 



3\(N) 



2 



^P-ldP-kaP-ka^ki+ka+ka + 7[7aA4 
ki,k 2 ,k 3 *' ' ki,...,k 4 



k 

1 

0,4 X 



P-k 1 P-k 2 P-k 3 P-k 4 ( 5ki+...+k4 , (17) 



where 



a n = + Aa n (18) 
and Aa n are the corrections obtained as the result of integration over CVs 
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c k : 



Aai = y^G(q), 



a a i 

A «3 = ^^EGfoWlq + k! |) + ..., 



fi S 1 

Aa 4 = ^_^G(g)G(|k 3 + k 4 -qi) + ..., (19) 



where 



^ = i ; \ (2°) 

is the charge-charge structure factor of the RPM determined in the Gaussian 
approximation. It is worth noting that G(q) is of the same form as the 
function G^^k) introduced in [TB]. In this study we restrict our attention to 
the GL critical point where G(q) remains a smooth function. 

Now let us consider the coefficient a-i- Expanding this coefficient at small 
k one can readily see that the linear term vanishes. As a result, we obtain 

«2 = ^2,0 + \k 2 a 2 ,2 + • • • , (21) 



2 



where 



1 ^G(a) d2d{q) (22) 
^A7)Z^)^2— ( 22 ) 



Inserting the expansion (|2Tj) into (fTTj) one arrives at the effective Hamiltonian 
of the RPM in the neighbourhood of the GL critical point. The coefficients 
of the Hamiltonian are given explicitly in (|14j) - (jl6j) and (|18j) - (j22j) . As is seen, 
the effective interaction appeared in (JT7j) has a short-range character. 
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4 Conclusions 



In this paper we have studied the critical behavior of the RPM at the vicinity 
of the GL critical point. For this purpose the method of CVs with a refer- 
ence system was used . First we obtain the functional of the grand partition 
function given in terms of the two CVs: the CV pk describing fluctuations 
of the total number density and the CV Ck describing fluctuations of the 
charge density. As is known, the Gaussian approximation of the functional 
of the grand partition function of simple fluids and their mixtures produces 
the qualitative picture of their phase behavior. In contrast, the full phase 
diagram of the RPM cannot be obtained within the framework of this ap- 
proximation. In order to describe the GL phase transition, the terms of 
the higher order should be taken into account in the effective Hamiltonian 
[THUS]. Actually, the charge-charge correlations causes the effective attrac- 
tion between the ions that, in turn, leads to the GL phase transition in the 
RPM. 

After the integration over CV Ck with the Gaussian basic density measure 
we construct the effective Hamiltonian in terms of CVs pk connected with the 
order parameter (up to p 4 ). All the coefficients of the effective Hamiltonian 
consist of two parts: the part depending solely on the characteristics of the 
RS (through the n-particle structure factors of the RS in long-wave length 
limit) and the part of the mixed type. The latter has the form of an expansion 
in terms of the charge-charge structure factors and is obtained as the result 
of integration over CVs Ck- Allowance for the charge-charge correlations 
(through the integration over CVs Ck) leads to the contribution Aa2 to the 
coefficient 02 (at the second power of CVs pk) which describes the effective 
attraction of the short-range character. 

Finally, the original Hamiltonian is mapped onto the GLW Hamiltonian 
of an Ising-like model in a magnetic field. We conclude that the form of the 
effective Hamiltonian of the RPM confirms the fact that the critical behavior 
of the RPM nearby the GL critical point belongs to the universal class of 
a 3D Ising model. The coefficients of the effective Hamiltonian are given in 
terms of the structure factors of the RS as well as in terms of the charge- 
charge structure factors. This demonstrates that the GL phase transition 
in the RPM is attributed to the effects of strong charge-charge correlations 
which take place at the low temperatures. A more comprehensive analysis 
of the coefficients will be done elsewhere. 
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Appendix 

We can rewrite (jlOj) in the form: 

J (p, c) = J (du) (drf exp ^i2n ^(w k Pk + 7kC k ) + 



2! 



J202 ] ^^ + 7k7-k) 



X 



1 + A + - { A 2 + 



where 



(23) 



n>3 i n >0 



We substitute into (J23j) the operators ^-^r and for cjk and 7k, re- 

spectively. As a result, we can write for J(p, c) 



J(p,c) 
where 



1 + i + -i 2 + . . . 



exp 



1 PkP-k 1 V-^ C kC-k 



- E 



E 



(24) 



i = EE^' l) ( 9 /^,^ Ck ). 

n>3 2n>0 

Substituting (J23j) into ©, we obtain for 5: 



S exp(/5pi07ti 0) ) y (dp)(dc) exp(/3pip 

lv-PkP-k 1 + 2$ C (A;)^ 2) 



1 + i + -i 2 + . . . 



x 



exp 



■ E 



2 V mtf 2V 



p y- 



anl 2) 



-CkC_k 



and after cumbersome algebra we can arrive at the expressions (I12j) - (jl6|) . 
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